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We use recent nuclear parton distributions, among them the Hirai - Kumano - Nagai 
(HKN) and Eskola-Paukkunen- Salgado (EPS08) parameterizations, in our pQCD- 
, improved parton model to calculate the nuclear modification factor, Raa'{pt), at RHIC 

and at the LHC. At RHIC, the deuteron-gold nuclear modification factor for pions, mea- 
sured at pt > 10 GeV/c in central collisions, appears to deviate more from unity than 
the model results. The slopes of the calculated RdAuipr) are similar to the slopes of the 
PHENIX pion and photon data. At LHC, without final-state effects we see a small en- 
^ l' hancement of RdPtiPr) in the transverse momentum range 10 GeV/c > pt > 100 GeV/c 

D I for most parameterizations. The inclusion of final-state energy loss will reduce the 

i-J^ . RdPbipr) values. 



I. INTRODUCTION 
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CO ' Nuclear modifications in the physics of partonic degrees of freedom arise either due to the 

ly-^ . environment in the initial nuclei around the incoming parton or from the interaction with the 

' partonic matter created in the final state of a nuclear collision. Coherent quantum effects may 

OO , connect and mix initial and final state phenomena. Here we investigate ^ cold nuclear effects^ 

' defined as effects that do not originate in connection with the thermalized partonic matter in 

collisions between two heavy nuclei [l|, [3| . 
• ■ Cold nuclear effects are usually described as modifications of the parton distribution functions 

K% ' (PDF) into nuclear parton distributions (nPDF). At sufficiently small momentum fraction, this 

^ . has been referred to as nuclear shadowing, but further effects, like nuclear multiple scattering, 

EMC, higher-twist effects, etc. need to be included in a more complete description. In this 
short review we present the high-p^ and high-xy behavior of some of these modifications, using 
shadowing/nuclear PDF parameterizations by HIJING d, EKS [J, EPS08 [|, and HKN d]. 

At present, cold nuclear modifications can be tested by experimental data at RHIC energies 
on hadron spectra from dAu collisions, or by direct photon spectra from dAu or AuAu collisions, 
which are essentially not influenced by non-Abelian jet energy loss. Note that recent theoretical 
studies point out the possibility of a small energy loss in pA or dA type collisions P, 0] . 
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II. COLD NUCLEAR EFFECTS AT RHIC 



We calculated the nuclear modifications at RHIC energies up topx ^ 70 GeV/c using various 
shadowing paramctcrizations in our pQCD improved parton model @] . Results are compared to 
the preliminary PHENIX data up to the highest measured momenta {pT ^ 18 GeV/c) [9, 10, 11]. 
On the left side of Fig. [TJ R^au^Pt) is plotted for neutral pions at -/s = 200 AGeV against 
a logarithmic pT scale, using the HIJING, EKS, EPS08, and HKN paramctcrizations. The 
HIJING shadowing is understood with its accompanying multiple scattering Q . 
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FIG. 1: Calculated RdAuiPT) compared to most central PHENIX data (left panel), and the extracted 



slope parameter, /? of Eq.([T} in different models (right panel) and Ref.s 



All calculated R^au^Pt) curves have an approximately linear dependence on log(pT) in most 
of the EMC region (15 GeV/c ^ Pt ^ 50 GeV/c). Thus, we used a simple linear function of 
log(pT) to get the slope parameter f3: 

RdAu{^Og{pT)) = « + log(pT) ■ (1) 

The slope (3 is plotted in the right panel of Fig. [1] We use the best-fitted linear curve in the above 
PT interval to compare the calculated slopes to the ones extracted from the data. The slopes 
of the model curves are quite similar, although they are smaller (in absolute value) than the 
measured quantities. This may be due to final state effects [l[ or isospin effects [l2| . In addition 
to pion production slopes in dAu collisions, direct photon production in both dAu and AuAu 
collisions is plotted. Our results are compared to the preliminary PHENIX data [l^, [Tl| . 
These data have decreasing tendencies with px, leading to negative P values similar to those in 
7r° production. In the case of the dAu photon data, the large error bars manifest themselves 
in a large uncertainty of the extracted slope. Direct 7 production in AuAu shows a stronger 
decrease, since the initial state nuclear effects enter the convolution integral twice in this case. 

We analyzed the theoretical uncertainties in our model [l|, 01 . We found that at pr > 3 
GeV/c our model is not sensitive to the scale choice. Nuclear PDF parameterizations vary 
both the position of the Cronin peak and the high-pT slope of RdAuiPr)- The HKN parameter- 
ization allowed us to analyze the errors via the Hessian method Q. An almost constant ±10% 
uncertainty was found in the whole transverse momentum range. 

It can be seen from Fig. [T] that the EPS08 and HKN nPDF parameterizations give slopes 
close to the experimental ones. However, HKN comes closer to the measured points at high pT 
(see the left panel) . Note that the agreement with the data can be improved by the introduction 
of a small opacity, as we have done in Ref. [l| . 



III. PREDICTIONS FOR THE LHC 



Based on Ref.s [l|, H, [l^ we estimated the cold nuclear modifications at LHC energies with 
larger intrinsic transverse momenta. We apply the same parameterizations used above. In 
Fig. [2]wc plot R^ji^{xt) in dPh collisions at midrapidity at RHIC and LHC energies as func- 
tions of transverse momentum fraction, xt = "ipr/y/s- The preliminary PHENIX dAu data 
on neutral pion production are also shown. The different panels display different shadowing 
parameterizations as noted on Fig. [2] 

Though at RHIC energies all parameterization give similar results, at LHC we see that at 
XT ~ 3- 10^"^ HIJING yields a strong (sa 30%) suppression, while HKN gives a 10% enhancement 
at the same xt values. The EKS and EPS08 results predict a 5—10% enhancement (suppression) 
around xt ~ 10~^ (10""^). The EPS08 results show somewhat stronger suppression than EKS 
at low XT, due to the inclusion of high-rapidity RHIC data in the fitting procedure of this recent 
parameterization. For the high momentum-fraction region both of these models show a similar 
behavior due to the EMC effect. HIJING and HKN have a reasonable agreement between RHIC 
and LHC energies (they scale with xt), while EKS and especially EPS08 display less precise 
scaling at around the anti-shadowing maximum. 




FIG. 2: Nuclear modifications -R^aCs^t) at RHIC (dashed) and LHC (solid) energies with different 
shadowing parameterizations. Preliminary PHENIX data are plotted as functions of xt- The effect of 
cold energy loss is indicated by thin dotted and dot-dashed Unes for L/\g — 1 and 3, respectively. 



Finally wc turn to possible cold energy loss at the LHC. The standard estimate is dN/dy w 
1500 — 4000 in central PhPh collisions at LHC energies [l3|- While in central dPh colhsions 
the geometrical cross section is small, using L/\ ~ {Npart^/^ from [| one obtains dN/dy « 
500 — 2000, still large, in central dPh collisions. This motivates us to examine the effect of a 
small energy loss in the final state. We apply opacity values L/\g < 3 for this illustration. The 
results are indicated as thin dotted {L/Xg = 1) and dot-dashed {L/Xg = 3) fines in Fig. [21 



The energy loss is stronger at the lower xt values. Here the relative particle yield is sup- 
pressed, and RdAixr) is steeper. In the high px region the suppression loses strength as it 
increases with ^ \og{E). We tested this effect at RHIC and LHC energies in AuAu (PbPb) 
collisions [THj . 

IV. SUMMARY 

We analyzed the high transverse-momentum behavior of the nuclear modifications in dA 
collisions. Using preliminary PHENIX data we checked several common nuclear shadowing 
parameterizations. We compared the logarithmic slopes of the nuclear modification factor in 
the 15 GeV/c Pt ^ 50 GeV/c region. The shadowing parameterizations investigated have an 
almost linear behavior with a negative slope in dAu collisions at RHIC, which can be attributed 
to the EMC effect. While we found all studied shadowing models reliable at intermediate pr, 
at high transverse momenta the EPS08 and the HKN nuclear PDF seem to give the best 
agreement with the PHENIX data. We tested the xt scaling of midrapidity 7r° production in 
dA collisions at RHIC and LHC energies in the models. At high xt most models have a similar 
behavior, but EPS08 deviates from precise scaling around a;^ sa 6 • 10^^. At lower xt, RdAixr) 
shows a 5 — 15% enhancement in the HKN, EKS, and EPS08 models. The enhancement is 
absent in the HIJING with multiscattering model. We tested the effect of a small energy loss, 
which counteracts any enhancement in dPb collisions and yields a wide-range suppression of 
the RdAi^r)- In light of these results and open questions, we emphasize the need for pPb (or 
at least dPb) measurements at the LHC. 
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